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Abstract 

Previous investigations into the morphology of Pyrrocoma liatriformis sensu la to in northern Idaho 
and adjacent Washington have revealed two distinct morphologies that correspond to their 
geographical ranges. These same populations and individuals have been analyzed using AFLP data. 

Over 400 loci were identified among all individuals using two sets of AFLP adaptors. The data are in 
agreement with the morphological data and separate the populations from the Snake River Canyon/ 
Camas Prairie from those of the Palouse grasslands. Data clustering methodologies using both 
presence/absence data for all individuals and allele frequencies for each population produced similar 
results. We suggest the name P. scaberula be resurrected to encompass the populations from the Snake 
River Canyon and Camas Prairie. 
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Numerous species concepts have been pro¬ 
posed to unambiguously determine species 
boundaries (see Niklas 1997; Howard and 
Berlocher 1998; Wilson 1999; Coyne and Orr 
2004; Sites and Marshall 2003, 2004, for more 
detailed summaries). Selecting a concept poses a 
challenge to biologists, particularly botanists, 
where relatively common gene flow among 
morphologically distinct populations seems to 
preclude the widespread adoption of the biolog¬ 
ical species concept (Burger 1975; Donoghue 
1985; Mishler and Brandon 1987; Ellstrand et al. 
1996). 

Over the past 20 yr, molecular methods have 
provided systematists an additional means of 
assessing species concepts beyond morphological 
variability (Miller and Spooner 1999; Lopez et al. 
1999; Duim et al. 2001; Parsons and Shaw 2001; 
Dawood et al. 2002; Wiens and Penkrot 2002; 
Richardson et al. 2003; Martinez-Ortega et al. 
2004; Sites and Marshall 2004; Whittall et al. 
2004; Garzon et al. 2005; Irwin et al. 2005; Pons 
et al. 2006; Suatoni et al. 2006; Roe and Sperling 
2007; Manoko et al. 2007; Meudt and Clarke 
2007; Guo et al. 2008). Molecular data allow a 
means to determine if populations are genetically 
distinct from each other, obtain estimates of gene 
flow between populations, and resolve if they are 
mutually monophyletic. 

Assessing population variation to resolve 
taxonomic status by molecular genetic means 


can be done by a variety of methods including 
simple sequence repeats (SSRs or microsatellites, 
Akkaya et al. 1995), inter-simple sequence repeats 
(ISSRs; Smith and Bateman 2002), randomly 
amplified polymorphic DNA (RAPDs; Smith 
and Pham 1996), amplified fragment length 
polymorphisms (AFLPs; Vos et al. 1995) and 
randomly amplified fingerprinting (RAFs; Wal¬ 
dron et al. 2002). Many studies have employed 
AFLPs to provide evidence that a single species 
should be divided into multiple species where 
morphological data were limiting or conflicting 
(Kardolus et al. 1998; Lopez et al. 1999; Mueller 
and Wolfenbarger 1999; Duim et al. 2001; 
Hedren et al. 2001; Koopman et al. 2001; Bottini 
et al. 2002; Parsons and Shaw 2001; Richardson 
et al. 2003; Martinez-Ortega et al. 2004; Whittall 
et al. 2004; Garzon et al. 2005; Irwin et al. 2005; 
Manoko et al. 2007; Meudt and Clarke 2007; Roe 
and Sperling 2007; Travis et al. 2008). According 
to Becker et al. (1995) AFLP fragments are 
derived from throughout the genome and bands 
of identical size (co-migrating) are predominantly 
homologous in closely related organisms (Waugh 
et al. 1997; Rademaker et al. 2000). Herein, we 
examine the AFLP variation among populations 
of Pyrrocoma liatriformis sensu lato. 

Pyrrocoma liatriformis E. Greene (syn. Hctplo- 
pappus liatriformis (Greene) St. John) is an 
herbaceous perennial found in northern Idaho 
and adjacent Washington (Fig. 1). This taxon has 


96 


MADRONO 


[Vol. 57 



Fig. 1. Map showing the distribution of populations sampled in this analysis. Abbreviations follow Table 1. 
Open circles represent P. scaberula and closed circles represent P. liatriformis, the star represents P. cartliamoides. 


generally been viewed with a broad species 
concept (Hitchcock et al. 1955) that may 
encompass two distinct morphological species. 
Recent work on the morphology of P. liatriformis 
sensu lato has discovered that neither qualitative 


nor quantitative morphological characters are 
uniform across P. liatriformis sensu lato (Bjork 
and Darrach 2009). Instead, the morphological 
variation is correlated with geographical range; 
plants from the Palouse grasslands are mostly 
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heavily tomentose throughout, resinous-punctate 
glands are absent and flower heads tend to be 
smaller (10-13.8 mm) whereas those of the 
canyon/Camas Prairie regions are hispid, with 
conspicuous resinous-punctate glands and flower 
heads are larger (13.8-15.1 mm). Taxonomic 
investigations of the type specimens of Pyrro- 
conia have revealed that these morphological 
entities have been described in the past as P. 
scaberula E. Greene (canyon/Camas Prairie 
plants) and P. liatriformis (Palouse grasslands); 
the former species largely has been considered a 
synonym of the latter in most treatments. 

Bjork and Darrach (2009) studied the mor¬ 
phological variation that separates P. liatriformis 
from P. scaberula by examining eight continu¬ 
ously variable morphological characters for 325 
plants, 201 individuals in 18 populations from the 
Palouse and 124 individuals in 13 populations 
from the canyon/Camas Prairie region. The 
characters showed clear non-uniformity and 
although there is overlap in the ranges of the 
characters, the means of six characters (lateral 
branches, number of heads, head length, head 
width, phyllary width, and leaf width) are 
statistically significant for each group of popula¬ 
tions. Furthermore, principle component analy¬ 
ses clearly separated the populations into two 
distinct groups. Although some populations in 
the Palouse grasslands share the lack of tomen- 
tum and strong glandularity of the canyon/ 
Camas Prairie plants, they were clearly attribut¬ 
able to the morphology of P. liatriformis sensu 
stricto based on quantitative characters. 

Despite the ability of the characters to separate 
the Palouse and canyon/Camas Prairie popula¬ 
tions into distinct clusters, there is a strong 
overlap among the ranges of the morphological 
characters. These data suggest the two species are 
either closely related, represent a progenitor- 
derivative species pair (Gottlieb 1973, 1974; 
Gottlieb and Pilz 1976; Crawford and Smith 
1982; Ranker and Schnabel 1986; Perron et al. 
2000), undergo hybridization, or are a combina¬ 
tion of these. 

It is the goal of this study to resolve whether 
the morphological species as defined by Bjork 
and Darrach (2009) are congruent with patterns 
shown by molecular data. 

Materials and Methods 

To obtain an estimate of molecular genetic 
variability, we sampled 32 populations of Pyrro- 
comci liatriformis sensu lato from both Palouse 
grassland and canyon/Camas Prairie populations 
(Fig. 1). One additional population of P. cartha- 
moicles was used as an outgroup for comparison. 
Populations, their abbreviations, and tentative 
species identifications are presented in Table 1 
and plotted onto a map in Fig. 1. These are the 


same populations and individuals sampled for 
morphological data by Bjork and Darrach 
(2009). At the time sampling was conducted, 
exceptionally dry conditions in the region result¬ 
ed in few populations that were flowering, which 
minimized the number of individuals that were 
sampled per population (Bjork personal obser¬ 
vation). Ideally 25-35 individuals were to have 
been sampled, but a total of 25 or more 
individuals was possible for only five of 33 
populations (Table 1). Leaves were collected on 
silica gel and were the source of DNA extraction 
using DNeasy kits (Qiagen, Valencia, CA). One 
fertile stem per plant from populations of larger 
size was collected as a voucher. Plants from 
smaller populations were vouchered nondestruc¬ 
tive^ with photographs taken with a ruler for 
scale. All vouchers, including photographic ones, 
are deposited at the University of Idaho Stillinger 
Herbarium (ID). 

Approximately 500 ng of DNA from each 
individual was digested with Mse I and EcoK\ 
while simultaneously annealing the Mse I and 
EcoPA adaptors at room temperature overnight. 
Annealing of the adaptors to the ends of the 
fragments alters the restriction site and precludes 
further digestion or need to perform reactions 
separately. To reduce the overall number of 
fragments amplified, and thus improve detection 
of homologous amplified fragments, a pre- 
selective amplification was performed using the 
primers £coRI + A (5'-GACTGCGTAC- 
CAATTCA-3') and Mse I + C (5'-GAT- 
G AGT CCT G AGT A AC-3') and 1 pL of digest¬ 
ed/ligated DNA. Pre-selective amplification was 
run with 20 cycles of denaturing at 94°C for 
30 sec, annealing at 56°C for 1 min and extension 
at 72°C for 1 min. Products were diluted 1:20 
with TE buffer and 3 pL of the dilution were used 
in selective amplification using either the A 
primer set, M-CAC (5'-GATGAGTCCTGAG- 
TAACAC-3') and labeled (with dye for Li-Cor 
system) E-ACT (5' -G ACT GCGT ACC A ATT - 
CACT-3') or the T primer set, M-CTC (GAT- 
GAGTCCTGAGTAACTC-3') and labeled E- 
ACC (5'-GACTGCGTACCAATTCACC-3'). 
Final selection amplification used an initial 
denaturation at 94°C for 2 min followed by 10 
cycles of denaturation at 94°C for 20 sec, 
annealing at 66°C for 30 sec and extension at 
72°C for 2 min. This was followed with 25 more 
cycles that differed only in reducing the annealing 
temperature to 56°C. Lastly there was a 30 min 
extension period. 

Final AFLP products were separated on 6.5% 
polyacrylamide gels and visualized on a Li-Cor 
LongreadIR automated sequencer (Li-Cor Bio¬ 
technology Division, Lincoln, Nebraska). Molec¬ 
ular weight size standards were run on each end 
of each gel. Digital images of the gels were 
analyzed using Gene ImagIR (Li-Cor Biotech- 
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Table 1. Locations and Abbreviations of Populations Sampled in This Analysis, Their Species 
Designation Based on Morphological Data, and Amplification Success for Each Set of AFLP 
Adaptors (A and T). Vouchers are deposited at the University of Idaho Stillinger Herbarium (ID). 


Population abbreviation Number of 

(number of individuals individuals amplified: 
Species designation sampled) (A/T) 


P. scabemla —Anatone, Asotin Co., WA 
P. scaberula —Chesley Railroad, Lewis Co., ID 
P. scaberula —Craig Mountain, Nez Perce Co., ID 
P. scaberula —Ferdinand Butte, Idaho Co., ID 
P. scaberula —Ferdinand east/Meadow Creek, Idaho Co., 
P. scaberula —Lawyer Canyon, Lewis Co., ID 
P. scaberula —Lime Hill, Asotin Co., WA 
P. scaberula —Redbird Road, Nez Perce Co., ID 
P. scaberula —Soldiers Meadow, Lewis Co., ID 
P. scaberula —Talmacks North, Lewis Co., ID 
P. scaberula —Upper Cold Spring Creek, Lewis Co. ID 
P. scabemla —Weissenfels Ridge, Asotin Co., WA 
P. liatriformis —American Ridge, Latah Co., ID 
P. liatriformis —Barking Dog, Whitman Co., WA 
P. liatriformis —Cedar Ridge, Latah Co., ID 
P. liatriformis —Eden Valley, Whitman Co., WA 
P. liatriformis —Genesee South, Nez Perce Co., ID 
P. liatriformis —Gross Road, Whitman Co., WA 
P. liatriformis —Joel, Latah Co., ID 
P. liatriformis —Kramer Prairie, Whitman Co., WA 
P. liatriformis —Lenville Road, Latah Co., ID 
P. liatriformis —Mix Road, Latah Co., ID 
P. liatriformis —Palmer Butte, Latah Co., ID 
P. liatriformis —South end of Paradise Ridge, Latah Co., 
P. liatriformis —Palouse Prairie Strip,Whitman Co. WA 
P. liatriformis —Rose Creek, Whitman Co., WA 
P. liatriformis —Armstrong Road, Whitman Co., WA 
P. liatriformis —Steptoe Butte, Whitman Co., WA 
P. liatriformis —Spaulding Road, Spokane Co., WA 
P. liatriformis —Uniontown, Whitman Co., WA 
P. liatriformis —Whelan Cemetery, Whitman Co., WA 
P. liatriformis —Wawawai Grade, Whitman Co., WA 
P. carthamoidesSmoot Hill, Whitman Co., WA 


nology Division, Lincoln, NE) to determine 
molecular weight designations for each fragment. 
Gel images were edited to ensure that fragments 
of identical size were correctly assigned the same 
weights. These data were exported using Gene 
Profiler (Scanalytics, Inc., Fairfax, VA) and 
fragments were assigned an allele designation 
based on which set of adaptors was used (herein 
designated as either the A or T set of alleles) and 
their molecular weight. Each individual was then 
scored for presence or absence of each allele. 

Fragments greater than 600 bp and less than 
49 bp were excluded from the AFLP analyses. 
Large fragments may be amplified with lower 
frequency during the process, as a result their 
consistency may be less reproducible and reliable. 
Similarly, homology of larger fragments becomes 
less probable since there are greater opportunities 
for insertions and deletions of DNA to alter 
fragment sizes between individuals. Smaller 
fragments were excluded due to potential diffi¬ 
culties in resolving fragment sizes. 


AN (12) 

12/3 

CH (25) 

21/24 

CM (24) 

24/23 

FB (25) 

22/24 

FE (16) 

13/13 

LC (12) 

11/12 

LH (25) 

22/25 

RR (25) 

23/25 

SO(l) 

1/1 

TA (12) 

0/12 

TS (5) 

5/5 

WR (8) 

7/7 

AR (16) 

16/15 

BD (16) 

16/16 

CR (10) 

10/10 

EV (7) 

6/7 

GE (4) 

3/3 

GR (11) 

10/10 

JO (5) 

5/5 

KS (19) 

19/18 

LR (11) 

10/11 

PP (5) 

5/4 

PB (5) 

5/5 

PR (15) 

15/15 

PS (9) 

8/8 

RC (5) 

5/5 

RT (6) 

3/2 

SB (7) 

111 

SP (3) 

212 

UN (3) 

2/2 

WC (25) 

20/23 

WW (21) 

20/20 

SM (23) 

23/22 


Data were entered into MacClade (Maddison 
and Maddison 2000) and exported as a simple 
table. The table was modified using AFLPDAT 
(Ehrich 2006) to convert the files to formats 
usable in STRUCTURE (Pritchard et al. 2000; 
Falush et al. 2007) as well as to generate allele 
frequency data for each population. 

There are many alternative approaches to 
analyze AFLP data to detect structure within 
the data set (Bonin et al. 2007). In general, these 
break down into two analyses: 1) analysis of 
bands directly (presence/absence) or 2) converting 
the band data into allele frequency data for each 
population. Both data types can then be used in 
an array of methodologies to detect diversity and 
structure within and among the sampled popula¬ 
tions. 

Here we opt to make use of both band data 
and allele frequency data to generate tree-based 
representations of the variation (Bonin et al. 
2007). We make use of two methods to analyze 
the data for both band and allele frequency data: 
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Neighbor-joining (NJ) using Jaccard’s genetic 
similarities (bands) or Nei’s (1972) genetic 
distance and UPGMA. These analyses are among 
the most widely used methodologies for resolving 
species boundaries using AFLP data (Miller and 
Spooner 1999; Duim et al. 2001; Koopman et al. 
2001; Parsons and Shaw 2001; Coulibaly et al. 
2003; Jacoby et al. 2003; Richardson et al. 2003; 
Dehmer and Hammer 2004; Martinez-Ortega et 
al. 2004; Whittall et al. 2004; Garzon et al. 2005; 
Manoko et al. 2007; Guo et al. 2008). 

A data matrix of presence/absence for all bands 
was directly imported into PAUP* (Swofford 
2002) for NJ and UPGMA analyses. Allele 
frequencies were used to calculate population 
genetic distances that were then used to generate 
NJ and UPGMA trees in PHYLIP (version 3.67; 
Felsenstein 2007). 

We used the Bayesian clustering method 
implemented in STRUCTURE 2.2 (Falush et 
al. 2007) to determine the optimal number of 
groups indicated by the data. Assuming all 
populations are in Hardy-Weinberg and linkage 
equilibria, this method assigns individuals to one 
of the pre-specified numbers of genetic clusters, 
K, using multi-locus genotypes and Markov 
Chain Monte Carlo sampling. We ran separate 
clustering simulations for all populations over a 
range of one to 40 clusters (i.e., K — one to 40). 
Individuals with data missing for either the A or 
T set of alleles were removed from the analyses. 
This reduced some populations to few samples 
(only three individuals for AN) and resulted in 
TA being removed completely. Simulations were 
run assuming an ancestry model that incorpo¬ 
rates admixture and correlated allele frequencies 
across loci for one million generations with a 
bum-in of 25,000 generations, sampling every 100 
generations. 

We compared posterior probabilities of K from 
one to 40 clusters using the ad hoc statistic, A K 
(Evanno et al. 2005). This statistic has been 
shown to be a better estimator of structure in 
some data sets, especially those where homoge¬ 
neous dispersal among populations cannot be 
assumed (Evanno et al. 2005; Travis et al. 2008). 
In such cases, a common pattern is for STRUC¬ 
TURE to plateau near the true value of K , and 
then to continue increasing gradually. Evanno et 
al. (2005) showed that A K consistently returns a 
clear peak at the true value of K under a variety 
of migration models. We chose a maximum K of 
40 because this exceeded the number of sampled 
populations. 

Results 

Amplification of DNA was successful for 
nearly all individuals for both sets of AFLP 
adaptors, A and T (Table 1). A few individuals 
did not amplify well which is perhaps the effect of 


plant resin that inhibits the reactions. No 
individuals from the TA population were ampli¬ 
fied for the A alleles and only three from AN 
were amplified fully for the T alleles. Some 
individuals were not amplified for either allele. 
Resin was detected in the precipitation stage of 
the DNA extraction for some individuals of these 
populations. 

For several of the gels, bands that were 
distinctly different in size below 50 bp were 
classified as the same size by the software. 
Therefore the limitations of the software for 
resolving bands at this stage required us to 
eliminate these bands from the analyses. Scan¬ 
ning gels visually indicated that these fragment 
sizes tended to be consistent across nearly all 
individuals, therefore their exclusion was unlikely 
to affect the results whereas their inclusion may 
have indicated erroneous relationships among 
individuals and populations. 

For the A set of adaptors, 245 alleles were 
scored for 373 individuals ranging from 569 to 
50 bp in size. For the T set of adaptors, 177 alleles 
were scored for 387 individuals ranging from 581 
to 50 bp in size. The complete data matrix had 
407 individuals from 33 populations and 422 
alleles. 

The analyses that used presence/absence data 
for all individuals were sensitive to missing data. 
Individuals with missing data for either allele 
were either 1) clustered together regardless of 
population designation, or 2) in disparate parts of 
the tree (often included in populations that were 
neither morphologically or geographically simi¬ 
lar, data not shown). Therefore we removed these 
samples (AN 1-9, AR8, CH2, 11, 12, CM11, EV5, 
FB9, 12, KS8, LC12, LH6,7,18, LR11, MIX1, 
RCA5, RR4, 7, RT4, TAOl-12, WC3,6,19,20, 
22,23,25, WW6, and 12) and re-ran the analyses. 
We also ran analyses of the A and T alleles 
separately to confirm the placement of the 
individuals excluded above in their respective 
population (data not shown). 

The analyses based on allele frequencies 
divided all populations into two distinct groups 
corresponding to Pyrrocoma scaberula and P. 
liatriformis based on a priori designations by 
Bjork and Darrach (2009; Fig. 2). Pyrrocoma 
carthamoides clustered within P. liatriformis with 
the UPGMA analysis (data not shown) because 
this analysis does not allow the outgroup to be 
specifically designated as such. The NJ tree in 
contrast (Fig. 2), results in two distinct and 
monophyletic groups each for P. liatriformis 
and P. scaberula. 

There are clear clusters within each of these 
groups. Populations that cluster together in both 
the NJ and UPGMA trees within P. scaberula are 
AN/FB/SO/TS, FE/WR, and LC/CH/LH/RR/ 
CM. Only population TA changes position 
between the analyses and is found close to CH/ 
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Fig. 2. Neighbor-joining based tree derived from AFLP allele frequency data. Population names follow 
abbreviations of Table 1. Bars to the right of the tree mark species boundaries, Psc —Pyrrocoma scaberula, Pli— 
Pyrrocoma liatriformis. Pea —Pyrrocoma carthamoides. Population names are abbreviated following Table 1 and 
are designated as normal font (P. scaberula ) or bold (P. liatriformis). 


CM/RR/LH/LC in the UPGMA analysis (data 
not shown). It should be noted that TA is lacking 
data for over half of the alleles. Within P. 
liatriformis there are similar clusters of popula¬ 
tions that are consistent between analyses. These 
are AR/PB/RC, PS/SB, JO/PP/PR, RT/WC, GR/ 
UN, and CR/GE/SP/KS/EV. The grouping of 
populations BD, LR and WW differ between the 
two analyses. With UPGMA, BD is close to all 
other populations of P. liatriformis, LR is close to 
RT/WC, and WW is close to CR/GE/SPKS/ 
EVLR/RT/WC (data not shown). 

The NJ analyses based on bands produced 
trees that were nearly identical to those based on 
allele frequencies (Figs. 2, 3) both in terms of 


clustering populations into species groups, and 
relationships of populations within each cluster. 
The greatest differences are that 1) not all 
populations were recovered as a single monophy- 
letic group (Fig. 3), 2) the NJ tree did not result 
in a monophyletic P. scaberula due to the 
position of population LH, and 3) P. cartha¬ 
moides is clustered within P. scaberula instead of 
P. liatriformis with the UPGMA tree (data not 
shown). 

Groupings within each species are also similar 
to the frequency-based methods. The major 
differences within P. scaberula are the position 
of population LH in the NJ tree (Fig. 3) and that 
individual AN5 was more closely associated to 







































































2010 ] 


SMITH ET AL.: PYRROCOMA SPECIES BOUNDARIES 


101 


— LH 


j—FE 
r P-WR 


-FE 12 


FE 13 


WR5&6 


AN5 


[C 


AN 


FB 


TS 


FB25 
-SO 


LC 


— RR 

-CM 

CH 


Psc 


r— PR 


Ljo 


pp 


SB 


PS 


.RC 

j- fLPB 


SB3&5 
-PS 9 


PS7 
PS8 


— AR 
-BD 


RC5 


TO 


WW21 


^LR 
-LR1&2 
L-WW7& 8 

f^GR 


UNI 


c 


c 


GR2,3 & 7 
GE 


UN2 


SP 


CR 


r WC 

-j—WC21 &24 

'-RT2 

ij—WC 4^,7^ & 18 

«-RT3 

EV 
—KS 


Pli 


SM 18 


I— SM 


Pea 


0.2 0.1 0.0 

Fig. 3. Neighbor-joining based tree derived from AFLP band presence/absence. Where the majority of sampled 
individuals formed a single cluster only the population name abbreviation is used. In instances where individuals 
fell outside of their respective population cluster, they are designated with a population name and number for the 
individual. Population names are abbreviated following Table 1 and are designated as normal font (P. scaberula) or 
bold (P. liatriformis). Bars to the right of the tree mark species boundaries, Psc —Pyrrocoma scaberula , Pli— 
Pyrrocoma liatriformis. Pea —Pyrrocoma carthamoides. 
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Fig. 4. Structure analysis of all populations of Pyrrocoma licitriformis sensu lato sampled in this analysis. A. Plot 
showing values of AK for each value of K. B. Bayesian assignment of individuals to two clusters. The bars represent 
the estimated posterior probabilities of each individual belonging to each of the two inferred clusters. 


WR/FE than to FB/SO/TS. The lack of unity of 
population AN is likely due to the fact that only 
approximately half of the T alleles were scorable 
for AN5. Population TA was excluded from these 
analyses due to a complete lack of A alleles. 
Based on analysis that utilized only the T alleles, 
TA was close to AN/FB/SO/TS as it is in the NJ 
tree based on frequencies (Fig. 2). Likewise, 
groupings were similar within P. licitriformis, the 
main exceptions being populations KS and EV 
which showed a closer affinity to WC/RT based 
on bands (Fig. 3). 

Plotting the actual values of K from one to 40 
indicated that there was a plateau after K = 2 
with small increases in probability for each 
subsequent value of K. With all populations 
except the outgroup P. carthamoides included, 
the greatest AK was at 2 distinct clusters (Fig. 4). 
These results agree with the clustering results that 
divide the populations into two species. Only 
individuals of population EV show any signifi¬ 
cant probability of being assigned to the other 
species (Fig. 4). 

Discussion 

All analyses of AFLP data presented here, 
regardless of whether bands or frequencies were 
used, separate the populations into Pyrrocoma 
liatriformis and P. scaberula as determined by 


Bjork and Darrach (2009) using morphological 
data (Figs. 2-4). The congruence of different 
methodologies is largely considered a means of 
overcoming potential problems of homology with 
AFLP data (Koopman et al. 2001) and the results 
of these analyses are congruent with previous 
work on morphology. 

Relationships Between Species 

Pyrrocoma scaberula is paraphyletic based on 
NJ analysis of bands (Fig. 3). This raises the 
question whether these two species may or may 
not represent a progenitor-derivative pair (Gott¬ 
lieb 1973, 1974; Gottlieb and Pilz 1976). The 
progenitor species would be expected to be 
paraphyletic since the derivative species would 
have resulted from a subset of populations. 
However, a second important criterion for a 
progenitor-derivative species is that the derivative 
species should contain a subset of the total 
diversity found in the progenitor. A summary 
of the presence/absence data shows that 73.2% of 
the alleles are shared between the two species, 
15.7% are unique to P. liatriformis and 11.1% are 
unique to P. scaberula. These results indicate that 
a large portion of the data is shared among the 
individuals and populations rather than being 
unique to either the putative progenitor species 
(P. scaberula ) or the putative derivative ( P . 
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liatriformis) and thus argues against a progenitor- 
derivative pair. Likewise, the results of the 
Bayesian simulations in STRUCTURE do not 
indicate any overlap of populations between the 
two species, but instead the optimal data 
partition is equivalent to two groups (Fig. 4). It 
seems more likely that the paraphyly of P. 
scaberula is a result of recent common ancestry 
with shared alleles between the populations. 

Relationships of Populations within Species 

The AFLP results clearly show population 
genetic structure within each species (Figs. 2, 3). 
Within P. scaberula, groups that consistently 
hold together following the tree-based methods 
include AN/FB/SO/TS/TA, FE/WR, and perhaps 
LC/LH/CH/RR/CM. Within P. liatriformis, AR/ 
PB/RC, SB/PS, JO/PP/PR, GR/UN, CR/GE/SP, 
and RT/WC are commonly recovered. Popula¬ 
tions BD, KS, WW, LR, and EV sometimes show 
relationships with other groups but not always. 
Population BD has the longest branches showing 
the greatest genetic distance from other popula¬ 
tions. The results from STRUCTURE indicate 
that the populations are best treated as two 
clusters, corresponding to the two species. The 
one exception is population EV where individuals 
show the greatest probability of being assigned to 
P. liatriformis, but also have a noticeable 
probability of being assigned to P. scaberula 
(Fig. 4). These data may indicate some recent 
hybridization within this population, be the 
retention of a large number of alleles that are 
ancestral to both species, or reflect incomplete 
lineage sorting. 
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